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Abstract 

A time-dependent approach is used to explore inelastic effects during electron transport through 
few-level systems. We study a tight-binding chain with one and two sites connected to vibrations. 
This simple but transparent model gives insight about inelastic effects, their meaning and the 
approximations currently used to treat them. Our time-dependent approach allows us to trace back 
the time sequence of vibrational excitation and electronic interference, the vibrationally introduced 
time delay and the electronic phase shift. We explore a full range of parameters going from weak to 
strong electron- vibration coupling, from tunneling to contact, from one- vibration description to the 
need of including all vibrations for a correct description of inelastic effects in transport. We explore 
the validity of single-site resonant models as well as its extension to more sites via molecular orbitals 
and the conditions under which multi-orbital, multi-vibrational descriptions cannot be simplified. 
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We explain the physical meaning of the spectral features in the second derivative of the electron 



g : current with respect to the bias voltage. This permits us to nuance the meaning of the energy 



value of dips and peaks. Finally, we show that finite-band effects lead to electron back-scattering 
off the molecular vibrations in the regime of high-conductance, although the drop in conductance 



I. INTRODUCTION 

The importance of inelastic effects in electronic transport in molecular junctions is widely 
recognized and it is a rich, active research field. Several recent reviews on the topic give 
a clear idea of its breadth^. Advanced new experimental techniques show that electrons 
transport through a few-atom system is strongly dependent on the vibrational degrees of 
freedom of the system. As miniaturization decreases device sizes, the role of atomic vibra- 
tions needs to be considered in the device functionalities. This dependence has been shown 
to profoundly alter the device behavior: from new channels of conduction^ to heating of the 
junction^. It is then important to understand the role of inelastic effects and the parameters 
that control them. 

A large body of research has been devoted to inducing controled inelastic effects. A recent 
review article^ shows that inelastic effects in a tunneling junction can be used to chemically 
analyze the molecules at the junction by inelastic electron tunneling spectroscopy (IETS) 6 , 
to induce reactions by displacing atoms 7 , and to use molecular conformational changes 
as a switch for possible devices^. In this way, the scanning tunneling microscope (STM) 
induces the reaction and also detects its product, a dramatic example is the modification and 
detection of trans- 2-butene on Pd (110)^. The effect of local vibrations during electron flow 
has also been revealed in careful experiments of photon analyses 10 yielding an unprecedented 
insight in vibrational dynamics. Not only has the tunneling junction of an STM been used 
to explore the coupled electron-vibration dynamics, also point contact spectroscopy has 
been used in molecular wires", where the conductance showed drops at the molecular 
vibrational onset. 

These many experimental results call for important theoretical development. Indeed, the 
last years have seen several theoretical works spanning most of the experimental systems: 
from tunneling inelastic spectra^^i to point contact spectroscopy^ 1 ^ 7 - 1 ^, from system- 
atic approaches studying different parameters' 1 ^ to lowest-order expansion with ab-initio 
parameters". A thorough review on methodology and results can be found in Ref. y. 

A quantitative description of inelastic processes is mandatory in order to assess the rel- 
evance of the different ingredients characterizing electron transport on the atomic scale. 
Recent developments in ab-initio calculations together with transport calculations permit 
us to grasp the essential parameters and, eventually, produce predictive calculations. Yet, 



the typical ab-initio-based calculations are on the one-hand-side heuristical, because a well 
established dynamical theory of electron transport is yet to come 23 , and on the other hand- 
side, they are complex and difficult to interpret. Most ab-initio approaches use ground state 
density functional theory with non-equilibrium Green's functions (NEGF), this combination 
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is not justified, and the results have the full complexity of NEGF. As signaled in Ref. 
new methodology to treat quantum transport will be probably based in time dependent 
density functional theory (TDDFT). Recent results show that it is possible to treat nuclear 
(semi classically) and electronic (quantally) degrees of freedom within TDDFT to treat the 
non-adiabatic transport problem 2 ^ 2 ^. New developments go a step further in the treatment 
of correlated electron-ionic dynamics^. But time-dependent methodology can have other 
benefits beyond its correctness. In particular, it can be used to develop a physical picture of 
the electron diffusion process, in this way yielding complementary information to the more 
involved NEGF approach. Time-dependent approaches can also have interesting numerical 
behavior. Indeed, electron transport treated with the short-iterative Lanczos method 2 ^ has 
a quasi-linear scaling for sparse Hamiltonians. 

In this article, we explore inelastic effects in electron transport by means of electronic 
wave packet propagations in an idealized atomic-size system. We consider tight-binding 
chains connected to one and two vibrating electronic sites. These vibrating sites can hold 
one and two nuclear modes that are coupled to linear order in the nuclear displacement with 
the electronic degrees of freedom. Despite the simplicity of this model system, the main 
one-electron ingredients are included, and the time resolved solution permits us to have 
insight different from the perturbation-theory Green's functions results. Similar treatments 
have already been performed for the case of electron-molecule collisions 2 ^ and for inelastic 
effects in transport 29 . 

The calculations presented here are distant from the experimental situation because the 
model system is very simplified and because many-electron problems are absent. Indeed, 
recent treatments show the richness of effects associated with the many-electron aspects 
of the problem 3 ^, as well as the non-equilibrium many-phonon problem 3 ^ 1 ^ 2 -. Despite, the 
absence of these very interesting ingredients of inelastic transport, our calculations can help 
in understanding inelastic effects because there are situations in which one-electron transport 
is justified even in the presence of inelastic effects^ 3 -. To a certain degree, our calculations 
are equivalent and complementary to those of H. Ness 3 ^ the main difference being that a 



time-dependent approach is adopted in the present study. 

II. TIME-DEPENDENT WAVE PACKET PROPAGATION 

Stationary electron transport does not need a time-dependent description. However, 
insight on vibrational excitation processes during the current flow can be gained by time- 
dependent calculations. Numerically, a time-dependent description can benefit from the 
quasi-linear scaling using sparse Hamiltonians. As in stationary descriptions, the bottleneck 
of the calculation lies in the matrix time vector product of the Hamiltonian acting on a 
system's vector. Hence, a time-dependent approach can yield especial insight in the actual 
transport process in a realistic nanoscale system by using a localized basis set that leads to 
a sparse Hamiltonian. This strategy seems to be particularly appealing for the description 
of conduction in nanowires^ and nanotubes^. In this section, we explore the general time- 
dependent methodology using wave packets, beyond the Kubo linear theory of Refs. I35IJ36 . 

A particularly interesting time-dependent approach to have information on the full elec- 
tronic trajectory is the short-iterative Lanczos (SIL) wave packet propagation^!. The initial- 
value problem is solved by applying the evolution operator to the initial wave function, so 
that the solution reads, *&(t) = e~* w '\l/(0), where Ti is the Hamiltonian of the system under 
consideration (atomic units are used throughout, h = m = e = 1, unless otherwise specified). 
The use of this equation is inconvenient in the case of large matrix dimensions because it 
implies a diagonalisation. Instead, we prefer to use a numerical approximation which con- 
sists in successive infinitesimal evolutions of the wave function, \I/(t + At) = e~ lAtH ^/ (t) . 
In order to have an efficient implementation, the SIL method truncates the Hilbert space 
to a subspace spanned by a few vectors. The computation of this truncated subspace is 
the bottleneck of the calculation because the matrix-time-vector product to generate the 
subspace is performed on the total dimension of the problem. Despite the truncation, the 
richness of the full Hamiltonian spectrum is recovered by repeating this operation for all the 
different time steps and propagating in this way the initial wave packet. 

Let $1 be a starting vector, the Lanczos propagation method consists in the construction 
of a Lanczos matrix following the recursion relation^ 1 ^, 

&+i$i+i = ^ $ i - «A' - &$J-i, J > 1 (!) 



where aij are the new diagonal matrix elements of the new tridiagonal matrix and (3j are 
the new upper and lower diagonal matrix elements. Since the recurrence relation is started 
by $1, initial wave function, \^(t = 0), then /?i = and <3?o = 0. The vectors <£>,, are called 
Krylov vectors, and define the Krylov space of order /, by spanning the subspace given by 
{$,-}, namely 

IC l = Span{^ 1 ,H^ u ...,n l - 1 ^ l }. (2) 

The Krylov vectors are orthogonal by construction, they may be normalized, and define 
thus a basis set in which the Hamiltonian can be expressed. In this basis, the Hamiltonian is 
tridiagonal. Additionally, the order / can be much lower than the initial matrix dimensions. 
Consequently, the Lanczos matrix can be easily diagonalized using conventional algorithms. 
Regarding the Krylov vectors, it is known that for relatively large values of /, their orthog- 
onality can be lost, this is the reason why we have set sufficiently small time steps to keep 
a rather small order I (I < 9). Its value can be changed at each step in order to optimize 
the performances of the calculation, lowering the computational cost as I decreases. An 
excellent description of this dynamical control of the accuracy of the Lanczos propagation 



method can be found in Ref. 
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Another important feature is the undesirable effects provoked by the finite size of the 
propagation grid. To eliminate artificial reflections of the wave packet at the boundaries, 
we use a parabolic absorbing potential. In order to account for the errors which may be 
introduced by the reflections at the boundaries in the presence of the absorbing potential, 
we calculated the reflection coefficient with a broad wave packet. We found that less than 
0.01% of the wave packet was reflected at any energy. 

To compute the energy-resolved transmission or reflection coefficients we use the virtual 
detector technique 40 which consists in the evaluation of the wave function some sites after 
the region where the electron- vibration interaction takes place. By a time-to-energy Fourier 
transform, transmission or reflection coefficients are obtained. In the inelastic case, this 
should be done for each vibrational state of the total wave function. Hence, for a 1-D 
tight-binding chain, the partial transmission, T n (oS), is given by 

nH ~ I « e M I 2 (3) 



where we have assumed zero temperature and an initial wave packet in the vibrational 
ground state, n — 0. ip^^uj) is the energy-resolved wave function in the vibrational state 
n, after the interacting region. The detector has been located on site d. ip^ a Q e (uj) is the 
wave function computed with the same virtual detector but considering a bare Hamiltonian 
containing no vibrational degrees of freedom nor elastic defects. 

To analyze quantities such as total transmissions, especially for non-trivial Hamiltonians 
which may have two or more vibrating electronic states, we have computed the density of 
states projected on any state, for instance, on any linear combination of the tight-binding 
states. Consider a state | a,n), with a a particular electronic state and n its vibrational 
state, evolving with Hamiltonian 7i. The density of states projected on | a, n) reads^*, 



00 dt 



p{uj)= —e lult (a,n\e- mt \a,n). (4) 

, Z7T 



III. THE PHYSICAL MODEL 

The Frohlich-Holstein model 42 is used to represent the combined electron-vibration sys- 
tem: the electron-vibration coupling is assumed to be linear in the normal-mode coordinates. 
If we assume only one mode of vibration, and a single site, the Hamiltonian reads, 

H = EqcIcq + Y^ £ k,ic{ yi C k ,i + Yl **,i(4,i C *+l.< + 4+l,i C M) (5) 

k,i k,i 

+ Y Mc^co + c l°k,i) + ^b + McJco(6 f + b) 

k.i 

where c k and Ck are the fermionic operators which create and anihilate an electron in state 
k. Similarly, W and b are the bosonic operators that create and anihilate a quantum of 
energy Q in the considered vibrational mode. The first term in the Hamiltonian ([6]) refers 
to the site where the interaction takes place, it has an on-site energy of e$. The second 
describes the energy of site k of chain i. For simplicity we will be dealing with only two 
chains (i = L,R: left and right). The third term describes the couplings among sites of 
chain i. Here, we just consider nearest neighbours. The fourth term is the coupling between 
the two semi- infinite chains and the state of energy eq. The fifth term is the energy of the 
harmonic oscillator. Finally, the last term describes the electron-phonon interaction. In 



this single-state single-mode model, M is only a scalar which represents the strength of the 
interaction. It is called the electron-phonon coupling. In the present work, we have explored 
both an electronic single site coupled to vibrations, and a double site. In the latter case, M 
will be a 2 x 2 matrix, as presented in section V. 

Using a tensorial product description of the electronic and nuclear coordinates, the full 
Hamiltonian can be expressed in matrix form: 
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where we have used a block representation in the vibrational basis {|n)}. The diagonal 
elements are electronic Hamiltonians, which define the propagation of a wave-packet in a 
vibrational subspace. For clarity we give their tight-binding representation, 
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where n labels a particular vibrational state in the harmonic approximation. The t are 
off-diagonal matrix elements which connect the nearest neighbours sites inside the left and 
right chain. The T\ end T r play the same role as t, they connect the chains to the site which 
has on-site energy eq. The diagonal term nfl accounts for the energy the electron must 
exchange with the vibrational degrees of freedom of the system. If it propagates from one 
vibrational state to another, it must lose or gain Q, the energy quantum of the vibration. 

The matrices M in (jSJ) couple the Hamiltonians H,^ in the different vibrational subspaces. 
In the case of a single-site impurity, M is essentially a sparse matrix with the same dimensions 



as H.^, where only one single element is non-zero, the one connecting the diagonal matrix 
elements of energy Eq + nVt and Eq + (n + 1)Q. In Hamiltonian ([6]), the factors that multiply 
M come from the matrix representation of the bosonic operators, they correspond to the 
factors which appear in the well known relations tf\n) = \/n + l|n+l) and b\n) = ^Jn\n — l). 
We note that Hamiltonian (^ is truncated, the number of vibrations that are considered in 
the calculation is N in this example. This number N may be sufficiently large to represent 
suitably the vibrational space. A detailed study of the convergency of the calculation is 
presented hereafter. 

Let us consider a molecule between two electrodes, modeled by a single state connected 
to two chains. In this single-site case, Y. Meir and N. S. Wingreen^ have shown that the 
current can be expressed as follows, 

J = -lJ[hH- f R (Lu)}Im{tr[TG r }}duj, (8) 

where T is defined as a function of the couplings to the right and left leads, V = f R +r ■ G r 
is the retarded green function of the molecule, and fuR) is the Fermi distribution of the left 
(right) lead. This formula applies in the case where the couplings to the leads only differ by 
a constant factor, A, such that I\ = Wr- This is always so in the wide-band limit for the 
single-site case. This is not the case beyond the wide-band limit. 

The equation above can be viewed as the integral of a quantity that we identify as a 
transmission, multiplied by an energy window given by the difference of the Fermi distribu- 
tions of the leads, meaning that a current will flow if both the transmission is non-zero and 
the voltage applied between the electrodes is sufficiently large. Following the derivation by 
H. Ness^ at the zero temperature limit, the transmission reads, 

Im{tr[T{uj)G r {uj)}} = r{uj)ImG r 00 , (9) 

where G r Q0 is the projection of the Green function in the n = vibrational subspace, Gq = 
(0|G r |0). This means that the transmission is related to the projected density of states, 
Eq. (jl|), in such a way that, if we consider the wide-band approximation, where the coupling T 
is independent of energy, the transmission is proportional to the density of states projected in 
the n = subspace. By using the optical theorem one can explicitly retrieve the vibrationally 
excited states in the inelastic current^. 



In the case of the wave packet calculation, we can extend these equations for the calcu- 
lation of the electronic current. We approximate the conductance assuming that the Fermi 
level only enters to define possible final electronic states (otherwise the calculation is fully 
one-electron, an in general we will not consider a Fermi level) 

<r(u;) = - J2 T n(^)lh(uJ + nQ) + f R (nQ - u)}, (10) 

n 

where the factor - is the conductance quantum in atomic units. The terms between brackets 
are introduced to account for the opening of vibrational channels since they contain the Fermi 
distribution functions of the left electrode, fi and the right one, fa, at zero bias voltage in 
the present case. Actually, when electrons do not have sufficient energy to deposit it into 
the vibrational degrees of freedom of the system, only the T (u) is to be considered in the 
conductance. The same holds if the energy of the incident electron is sufficient to excite one 
vibration, in this case we consider the sum of %(uj) and T\{uj). Nevertheless, the calculation 
should be performed with a sufficiently high value of n in order to take into account the 
influence of closed channels in the conductance. Finally, the current as a function of the 
voltage, V, can be written as the integral over energy of the conductance, 

/oo 
<r{u)[f L {oj)-f R {u)]dw. (11) 

•oo 

The terms between brackets ensures that the transmitted electrons go from occupied states 
to empty ones where the voltage dependence is included. 

The treatment presented here is complementary of the one by H. Ness^. Indeed, the 
physical model is the same one, the difference is the solution method. As in the case of 
Ref. [34| the present approach is single-electron, neglecting both electron occupation and 
electron-electron correlation effects. 

IV. SINGLE-SITE RESONANT MODEL RESULTS 

Let us assume a single state coupled to the two 1-D tight-binding chains of the above 
model. Figure Q] shows the results of a wave packet propagation. Figure Q] (a) is the elastic 
wave packet, n = 0, because we assume that no vibration is initially excited in the system and 
the temperature is zero. The present system is two 1-D electodes symmetrically coupled to 
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FIG. 1: (Color online) (a) Spatial dependence for the squared modulus of the electronic wave 
function at different times, arbitrarily shifted for representation purposes. Case (a) is the elastic 
propagation of the wave packet, n = 0, where n is the vibrational quantum number. Symmetric 
hopping matrix elements, Tj_, = Tr have been used. In (b) the wave packet is represented at the 
same instant in time for different vibrational states n. The wave packet is initially propagated from 
the left electrode in the vibrational ground state, n = 0. The calculations are performed at zero 
temperature. When the wave packet starts populating the site number 600 (marked with a cross in 
(a)), the different vibrational states are populated (b), and the vibrationally excited wave packet 
starts propagating in both directions, indistinctly, due to the symmetric electronic couplings. 

the vibrating site. Hence, in Fig. CD (b) we see that only the reflected wave packet is different 
from the transmitted one in the elastic channel, n = 0, while it is identical for the inelastic 
channels: the electron reaches the active site and populates the ladder of excited vibrational 
states. In each state the electron has a finite and identical left and right transmission. 

A simple-minded picture of the structure appearing in the transmission as a function of 
incident electron energy could be that there is a series of N resonances displaced by the 
phonon energy, in agreement with the scheme of Fig. [2j At a given inital electron energy, 
the wave packet probes the resonant electron site if the energy is within nfl of the resonant 
site energy. The results is a series of equidistant peaks spaced by Q. This is plotted in 
Fig. [3] (a). There, the transmission for the single electronic state without electron- vibration 
coupling is depicted in a dotted line, and the full line corresponds to the case when the 
electron-vibration coupling is included. 

However, the physical picture is more complex than just a series of equidistant reso- 
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nances. In order to understand the appearance of the vibrational sidebands we divide the 
transmission according to the final state of the vibrator, Fig. [3] (b). This yields information 
on the vibrational state once the wave packet has propagated through the resonant site. 
The result is that each transmission curve for a well-defined final vibrator state displays 
a similarly rich peak structure. The above picture, Fig. [2, has to be changed: there are 
complex vibrational pathways, in which different vibrational states are probed before the 
system is left in a singly well-defined vibrational state. Given the coherence of the electron 
propagation, the different pathways can interfere and will give rise to Fano lineshapes in the 
transmission function^. 

A. Vibrational state population sequence: electron coherence 

Wavepacket dynamics can probe the different structures appearing in the transmission 
function in order to yield information on the actual interference process. Let us take a 
vibrational ground state, broad electronic wave packet energetically centered about the first 
maximum of the full transmission function, Fig. [HJ A spatially broad wave packet is almost 
a planewave and synonymous of a monochromatic packet, hence we can be sure to probe 
only the structure of the first maximum. 

Figure H] shows the modulus square of the electronic wave function for the first three 
vibrational levels on the vibrating site as a function of time. We first notice that the 
population of the vibrations are sequential: the peaks are slightly shifted in time as n 
increases. The linear electron-vibration coupling forces this sequential population since n 
is changed in steps of one quantum of vibration, Eq. ©. The populations show a time 
modulation, in the present case, we see that the modulation frequency of the n = 1 wave 
packet is roughly twice the modulation of the ground state and of n = 2. This is easily 
explained by considering populating and depopulating the n = 1 level by populating the 
n = and the n = 2. After a certain time the n = 2 level can depopulate again in the 
n = 1 as well as the n = 0. Despite the simplicity of the electron- vibration coupling and 
the electronic model, the final population dynamics depend on the wave packet energy, the 
strength of the vibrational coupling and the electron lifetime at the site, leading to non- 
trivial dynamics. The lifetime of the site is fixed by the hopping matrix elements 7} and T r , 
Eq. (ffl). 
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FIG. 2: Vibrational energy level scheme: The electronic wave packet propagates from left to right 
in the n = vibrational ground state. When the electron wave packet populates the impurity site 
connected to the vibrator, the population of n = l,n = 2, ..., becomes different from zero and the 
wave packet propagates in both directions. The transmitted wave packet permits us to compute 
the transmission resolved in n. The transmitted wave packet can be assigned to a vibrational 
channel. The n = channel is the elastic one, that can have contributions from all channels due to 
the excitation and de-excitation of the vibration during the wave packet propagation. This leads 
to a rich vibronic structure in the electron transmission even for the elastic channel. 

The Fano profile characterizing the transmission functions, Fig. [3l can then be explained 
by the interference of several of these vibrational pathways. Indeed, asymmetric Fano profiles 
are the rule in Fig. [S] (b). For the case of the total transmission, the addition over final 
vibrational states smear out the different profiles, rendering more difficult the determination 
of an asymmetric Fano profile in the peak sequence. 
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FIG. 3: (a) Electron transmission for two 1-D tight-binding electrodes connected by a single 
electronic site with electron-vibration coupling (full line) and in its absence (dashed line), (b) 
Tranmission decomposition in the final vibrational state of the single site, n is the number of 
vibrational quanta. The elastic transmission is the n = curve in (b). The initial state Eq = —0.7 
eV, ft£l = 0.05 eV, the electron-vibration coupling is M = 0.04 eV and the resonance width is set 
by the hopping matrix elements T[ = T r = 0.05 eV. The bandwidth is 2 eV. 

B. Time delay and phase shifts 



Spatially narrow wave packets contain most of the energy components needed to span 
the spectral region of interest, i. e. the main peaks of Fig. [3J The Fourier coefficients in 
the energy domain can be evaluated to obtain the phase of each component and hence the 
phase shift between Fourier components. The derivative of the phase shift with respect to 
energy yields the time delay of the Fourier component in the wave packet^. Hence, we can 
analyze the effect of the vibrational excitation of the electronic site by means of the time 
delay imposed on the propagating wave packet. 

Figure 0(a) shows the phase shift for the elastic component of the wave packet (the n = 
component). The full phase shift over all of the vibronic peaks is n, this is the phase shift of 
the electronic resonance. The phase shift of each individual contribution to the transmission 
is more complex. We see that the phase shift varies rapidly after each maximum, very 
much indicating a suite of Breit-Wigner-like resonances. Indeed, the the phase shift 5 for a 
Breit-Wigner resonance of FWHM T centered at E is given by~: 

r 



5 = —tan 1 



2{E - E 



13 



0.3 

0.25 

0.2 

£10.15 

0.1 

0.05 





1 ' ' 


1 




1 


, | i , 


1 ' 


Va 






(a) : 


1 


i \ 

\» \ 

\ \ \ 






-11=0 
-n=l 


~ ■' 


\ M 






- n=2 


III 


1 ' 






■ 




'wV.\ 

1 


1** 


*^N-*^ __ 







2000 



4000 
Time 



6000 



8000 



0.3 

0.25 

0.2 
n 

010.15 

0.1 

0.05 





1 ' ' 1 


1 ' i 


ll 1 1 , 


A 




(b) . 






' 


I A w\ 




-n=0 


r A 




-n=l 


- .'i \ i 




- n=2 


■ ■•' W ' 




■ 


- »i \/i>. 


( X A 


- 


" v i 




. 


- fe V 


'VA \ 


- 


/,/ •' 


\ .\,A\ 




, Jf , , 


7 •V\-» 

I'- 


C"»^N. 







2000 4000 6000 8000 

Time 



FIG. 4: Time dependence for the squared modulus of the electronic wave function at the site 
connected to the vibrator. In case (a) an almost monochromatic wave packet has been centered at 
the energy of the lowest peak of Fig. [3] (b). In the case (b) the wave packet has been centered at 
the second peak, where the transmission in n = 1 is smaller than the transmission in n = 2. 

The time delay, r, is given by the energy derivative 28 of the phase 5: 

_ d5_ 
T ~dE' 

In the case of a single Breit-Wigner resonance, the time delay is then 

r/2 
T ~ (E-E y + (r/2y 

On resonance the time delay is just r = 2/T, yielding direct information on the resonance 
width, T. In the present case, the time delay, Fig. 0(b), gives some definite interpretation: 
the time delay is maximum once the wave packet has encompassed one of the vibronic 
resonances. Hence, just above the resonance, the electron is deterred in its propagation 
by interaction with the vibration. The vibration contributes to the partial width of the 
electronic resonance, however the total width is independent of the vibration 47 . A spatially 
narrow wave packet, will have the phase shift and time delay of an electronic resonance 
regardless of the existence or not of the vibrations. We also find negative time-delays, that 
correspond to drops in the phase shift. Electrons can hence be expelled from the resonance 
more easily in the presence of vibrations at certain incident energies. This behavior is due 



to the interference between two consecutive vibronic resonances. As emphasized in Ref. 



23, 



the time delay has to be interpreted with care in the case of a multichannel problem such 
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FIG. 5: Electronic phase shift (a) and time delay (b) for the elastic channel, n = 0, for the case 
of Fig. [3] (b). The phase shift over the full resonant structure amounts to tc, but the vibrational 
substructure leads to fast variations of the phase, showing the resonant nature of the vibronic 
peaks of Fig. El The time delay is maximum away from each resonance, but it can change signs on 
resonance, depending on the order of the vibronic peak. This delay/acceleration process is related 
to the interference path inside the vibrator. 

as the present one. The time dependence is indeed complex and the time delay is just a 
number that cannot summarized the full intereference pattern. 

The interferences among vibrational paths are ubiquituous in all present results. Their 
effects can be seen in the Fano-like lineshapes of the tranmission function Fig. [31 in the 
oscillating population of vibrational states with time, Fig. 0] and in the rapid changes of the 
phase shifts and in the negative time delays, Fig. [51 

C. Time scales 

To develop an understanding of the vibrational excitation process, many authors resort to 
comparing the different timescales in play^ 1 ^. This is a very appealing approach because 
it permits us to rationalize the excitation process and the excitation regime of different 
systems. 

Let us define T mo i as the lifetime of an electron in the molecule, which is given by the 
inverse of the resonance width, 1/T. Typically, for chemisorbed molecules, T ~ 1 — 3 eV, 
which means electron lifetimes in the sub-femtosecond range. In order to estimate the 
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change in conductance due to an inelastic process, we can compute the inelastic fraction 
of electrons passing by the molecule. This can be estimated by computing the ratio of the 
excitation time to the first quantum of vibration, T n= x, and the lifetime of the electron in 
the molecule, 7" mo i. Let us estimate the magnitud of the electron-vibration coupling, M, 
to attain a measurable change in conductance (larger than 1%). This is the regime where 
perturbation theory is valid. In this case, the excitation process has been classified as sudden 
in the electron-molecule collisions literature^. 

In a sudden process, the molecule is assumed to be very briefly in its negative ion potential 
energy surface (PES). This can be indeed very brief as has been said above. As in Ref. y, 
let us assume that the PES is basically parabolic, but displaced with respect to the neutral 
PES-^: 

E_= l -K{Q-Q_f (12) 

where K is the spring constant of the PES related to the mode frequency by K = fiQ 2 
where /x is some reduced mass (this is easily generalized to the case of multinuclear modes). 
Here, Q- is the displaced center of the negative PES. When the electron flows through the 
molecule, the molecule is suddenly in its negative PES. Hence, the nuclei experience a force 
given by 

F- dE ~\ 



Then, they acquire a speed, v, of the order of 

ldE_ 



and T mo j « 1/r. Hence, 



H dQ 
KQ 



q=o x r mol 



/xr ' 

where KQ_ is the electron-vibration linear coupling, M/5Q rms of Eq. (jBJ) with 5Q r 



l/\/2/ifi coming from the second quantization of the displacement Q. The speed gained by 
the nuclei at an excitation of one quantum of vibration near Q = is vm = y2f2//L In this 
way, we find an upper limit for T: 

r<|M|. (13) 

This simple estimation shows that for strong electron-vibration coupling, vibrational exci- 
tation is unavoidable, where strong means larger than the molecule-electrode coupling. 
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Weak coupling is then the regime when \M\ <C I\ We can use Fermi's golden rule to 
estimate the vibrational excitation rate 55 : 



T n =l r 

Assuming typical IETS branching ratios 7Vriol/ T ii=l ~ 10~ 2 , leads to 

\Ml < J_ 

r 2 - 400' 

therefore the coupling becomes \M\ ~ T/200. For chemisorbed molecules this is in the 
range of 0.01 eV. This coupling is easily found in a large class of molecules, and therefore 
IETS is a feasible spectroscopy This is a simple estimation showing that rather small 
M (~ 0.01 eV) can have a measurable change in conductance for molecules adsorbed on 
metallic leads. 



One more time scale is fixed by the vibrational frequency, Q. As discussed in Ref. |28 , 
there are three basic regimes that can be distinguished: 

1. r <C fi the negative ion is long lived and the molecule can vibrate in the new state, 

2. T 3> Q the negative ion is short lived and we are in the above sudden regime, 

3. T m Q strong interference effects appear due to the nuclear motion. 

The time scale given by Q, defines then the type of electron-vibration interaction that will 
result. In the first case, the electron has ample time to interact with all the nuclear degrees 
of freedom and depending on the electron-vibration coupling, M, vibronic signatures can 
appear. This has been shown in the case of STM studies of electronic states on surfaces^ 5 ^ 
with a band gap leading to small T. The second case corresponds to the IETS case discussed 
above. The third case has been studied in the gas phase (see for example Ref. 1281 ) but we 
are not aware of any report on the consequent boomerang effect in the electron transport 
regime. 

D. Finite-band effects 

The wide-band approximation simplifies the expression for the transmission that becomes 
analytical 47 . One of the results of the wide-band approximation is that the transmission 
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function and the spectral function or PDOS, Eq. (j3j), have the same behavior with the 
electron energy. This is defmetly not the case when finite-band effects are considered, indeed, 
the energy dependence of the coupling to the electrodes leads the transmission function to 
have a different energy behavior than the site spectral function, Eq. ([9]). When the initial 
level, Eq is near the center of the band, the conditions for the wide-band limit are easily 
attained. As the level approaches one of the band edges, T becomes strongly dependent of 
energy as well as the center of the transmission distribution, Fig. [31 

One of the apparent features is that the transmission peaks become thinner and better 
defined. As Eq approaches the band edge, many inelastic channels start closing because 
the final electron energy falls out of the electron band, leading to a substantial decrease in 
the transmission peaks. This is clearly seen in Fig. EJ The parameters of the transmission 
calculated in Fig. [U]are exactly the same as for Fig. [3J (a) except that e is -0.95 eV instead 
of -0.7 eV. The bottom of the band is at -1 eV and the quantum of vibration is Ml = 0.05 
eV. The singular aspect of the peaks is enhanced by the 1-D character of the electrodes. 
Indeed, when the final energy of the electron approaches the bottom of the band a van Hove 
singularity appears in the density of states making the transmission more singular. In more 
realistic cases with 3-D electrodes we expect less singular transmission peaks. 

For an electron resonance approaching the top of the band, the situation is rather different. 
Here, the higher vibronic side bands dissapear but the n convergency remains the same 
because no channel is closed, leading to low-energy side bands of the same width as for the 
case of the resonance at the center of the band. 

E. The dynamical polaron shift 

The polaron shift is defined as the energy displacement of the first peak in Fig. [3J (a) 
with respect to the dashed-line peak. The polaron shift is related to the appearance of a 
vibronic structure. As we showed above, in the time-scale discussion, 2 conditions need be 
satisfied. Namely, M must be sizeable and r <C Q. This leads to considering^ the parameter 
p = M 2 /(TQ). When p > 1 the polaron shift becomes measurable. 

Hyldgaard et al.— show that the polaron shift critically depends on the initial occupation 
of the vibrating electron site. When the electron site is occupied, the transmission function 
shows a main peak displaced by — 3M 2 /Q. At half-filling the displacement is —2M 2 /Q and 
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FIG. 6: Transmission function for the same system of Fig. [3] where Sq is -0.95 eV instead of -0.7 
eV. The closing of channels due to the band edge leads to sharper and singular-like transmission 
peaks. The band edge is at -1.0 eV. 

for an empty site the polaron shift is —M 2 /Q. 

We can give a direct interpretation of the polaron shift by using the above parabolic 
model, Eq. (1121) . Following the above discussion, the polaron shift for the empty site is 
—M 2 /Q which is equal to — \KQ 2 _. The neutral PES is given by ^KQ 2 , hence the polaron 
shift is exactly the energy difference between the neutral PES and the negative one at the 
nuclear coordinate at which electron capture takes place. Hence, the polaron shift is the 
energy gain in the formation of the negative intermediate because the electron-vibration 
coupling permits the electron to probe the negative ion PES. In the absence of coupling, the 
nuclear coordinates do not evolve, and the negative ion resonance is just a simple lorentzian. 

In order to probe the polaron peak, the nuclear wavefunctions in the two ground states 
(neutral and negative) need to overlap. The overlap is large at weak electron-vibration 
coupling, hence a well-defined peak appears, just shifted by the energy gain. As the coupling 
increases, the overlap diminishes, leading to a decrease of the polaron peak. In the limit of 
strong coupling, the polaron peak is basically zero and a rich structure of evenly distributed 
peaks of the negative-ion vibrational states is apparent, with a gaussian envelop^ 1 ^. 
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FIG. 7: (Color online) Convergency in the polaron shift and inter-peak distance with the number 
of vibrational states included in the calculation. The polaron shift is the displacement of the first 
peak with respect to the resonance without electron- vibration coupling, dashed lines in Fig. [3] (a). 
It is a critical measurement of the strength of the electron-vibration coupling, and hence, sensitive 
to convergency. Also sensitive is the inter-peak distance that is equal to S7. 

The polaron shift is a good test for the convergency of numerical calculations with the 
number of vibrational states included in the expansion of Eq. ([6]). Figure [7J shows the 
convergency of the polaron shift with the number N of phonons included in the calculation. 
A correct value is obtained (~ — 0.038eV) at about N=5 phonons. Following the results of 
the calculation of Hyldgaard et al. for an empty state, the polaron shift would be — 0.032ey 
with our parameters (M = 0.04eV, Q = O.OSeV) 

We also show in Fig. [7J the inter-peak distance as a relevant quantity regarding the 
convergency of the calculation. We observe that the peaks of lower energy tend to converge 
with less phonons than the peaks of higher energy, which need a higher number of phonons 
to be placed at Q the one with respect to the other. Trivially, higher-order peaks need higher 
n, otherwise the peak may not even appear in the calculation. 

The polaron shift is also very sensitive to the different approximations that are used 
when evaluating inelastic effects in electron transport. Indeed, the self-consistent Born 
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approximation (SCBA) 3 — yields wrong polaron shifts. Such a study is perfomed in |34. 
There, it is shown that the SCBA is equivalent to neglecting the ^/n factors of Eq. (jSJ), this 
leads to wrong interpeak distances. Nevertheless, the SCBA captures much of the physics of 
inelastic processes and is an all-order theory, becoming reliable enough and very interesting 
for the evaluation of inelastic processes in a wide range of problems. In the particular case 
of weak electron- vibration coupling, n « 1 and the SCBA is virtually exact. 

V. TWO-SITE RESONANT MODEL RESULTS 

Most of the literature devoted to transport in the presence of electron- vibration coupling 
is based on the single-site model. However, Bringer et al.— showed that the single-site case 
has a behavior that cannot be extended to more realistic systems in which several electronic 
states are coupled with vibrations. Here, we will analyze when one can reduce the problem 
to the single-site case and when not. In the same way, we will show that when several 
vibrations are involved, the neglection of some of the vibrations can lead to qualitatively 
wrong results. 

A. Two-sites and one vibration 



Our previous one-site model is extended to have two electronic sites connected to a vibra- 
tion. This model has been recently explored in the context of ab-initio based calculations 
of inelastic transport between two pyramidal electrodes 5 - 9 -. The simplified two-site model 
permits us to understand the more complex ab-initio results. The model is given again by 
Hamiltonian ([6]) and (EJ) where Eq and M are 2x2 matrices. 

Here, one only vibration is assumed to interact with the electron flow. In the case of weak 
coupling 59 , this assumption is justified because Eq. (JHJ) is truncated to the first phonon of 
each mode and the inelastic effects of each mode are separable. In the next section, we will 
see that when convergency in the number of phonons needs more than one excited vibration, 
all of the vibrations need to be considered at once. 

In this one-dimensional model, there are two only possible modes^ 9 -. We will call them the 
center-of-mass mode (CM), and the anti-bond- length mode (ABL). The CM mode means 
that both sites displace in phase, hence the electron- vibration coupling matrix have non-zero 
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FIG. 8: (Color online) Electron transmission as a function of incident energy for the two-site 
problem in semilog scale. Two different modes are consider, the anti-bond-length modes (ABL) in 
full line with dots, and the center-of-mass mode (CM) in dashed lines. These two calculations are 
performed independently, however when both modes are simultaneously considered the transmis- 
sion function is different and becomes the full line marked by ABL & CM modes. We see that no 
direct assignments of the transmission peaks can be performed and the vibronic structure grows 
in complexity beyond a superposition of transmission peaks from both modes. 



onsite elements (the diagonal): 



M 



CM 



m 3 
— 777.3 



(14) 



and the second element is negative in order to account for the sign of the displacement of 
each site. In the ABL mode, the sites are moving out of phase, corresponding to an internal 



stretch mode. Hence, the electron-vibration matrix becomes 



M ABL 



mi rri2 
m 2 mi 



(15) 



Figure [H] shows the transmission for the ABL and CM modes, independently computed 
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FIG. 9: Projected density of states (PDOS) on molecular orbitals: bonding (<r) in full line and 
antibonding (a*) in dashed lines. The modes are considered separately in order to analyze the 
transmission functions of the ABL mode (full line with dots in Fig. [8]) and of the CM mode 
(dashed line in Fig. [8]). In (a) the vibronic density of states for the ABL mode is projected onto 
the two molecular orbitals and in (b) the vibronic density of states of the CM mode. In (a) we see 
two distinct series of vibronic peaks corresponding to the bonding and the antibonding molecular 
orbitals since the ABL mode does not couple these molecular orbitals. However (b) shows mixed 
structure due to the coupling of molecular orbitals by the CM mode. 

and together (to be analyzed in the next section). These figures show series of vibronic 
peaks difficult to analyze. 

In order to analyze the electron transmission, Fig. [8], it is more convenient to use projected 
density of states (PDOS) of the full vibronic structure on molecular orbitals, Fig. [5](a) for 
the ABL mode and (b) CM mode. These orbitals are linear combinations of the two sites: 



and 



\°) 



IO 



^\L) + \R)) 



^(\L)-\R)) 



that diagonalize the uncoupled two-site Hamiltonian with eigenvalues e — £ mo i and ^o+^ mo L 
respectively, where Eq is the level energy for each site and t mo i is the hopping matrix element 
between sites. 
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In this new basis set, the above coupling matrices, Eq. (fT4"|) and (IT5|) . become: 



1713 
m 3 



(16) 



and 



M 



ABL 



(17) 



mi + rri2 

mi — rri2 

that hint at the interpretation of the above vibronic peaks. In the case of the CM mode, the 
a* and a orbitals are coupled via the vibration, while in the ABL case, the molecular orbitals 
are not mixed by the electron-vibration interaction. This case can then be interpreted as two 
single-sites connected to a vibration, and hence two vibronic sequences are associated with 
the spectral feature of a* at and the spectral feature of a at . Since, the effective electron 
vibration coupling is mi + rri2 for a, the number of peaks and the general vibronic sequence 
corresponds to stronger coupling than the vibronic sequence of the a* peak. This is clearly 
seen in the PDOS on a* and a, Fig. [HI where we see that the vibronic sequences with a* 
and o character are energetically localized near the original molecular-orbital peaks. 

The CM mode mixes a* and a, hence we obtain vibronic peaks that are shared in the 
PDOS over both molecular orbitals. 

For small enough couplings, we can estimate the vibronic structure keeping just n = 1 
in the vibrational part. Hence we can diagonalize Hamiltonian ([6]) for the sites connected 
to the vibrations. This can also be expressed by a hybridization scheme, Fig. ffTUl) . We 
realize that the new peaks correspond one-to-one to the peaks found in the PDOS, and that 
they have contributions in different ratio from the a* and a orbitals according to the matrix 
elements of M. Depending on the magnitude of qm = ifr 2- ) 2 more or less phonons, n, will 
be needed to converge the vibronic sequence^, and more or less elements will be included 
in scheme f JTQ|) . In order to obtain the type of diagram (a) or (b), only knowledge of the 
symmetry of the system is required. 

In the present case, while the CM mode mixes a* and a orbitals efficiently, the ABL mode 
does not mix them. Hence, only this last case can be understood by a single-level model. 

It is interesting to study the time dependence of the electron occupation of both sites 
depending on the considered mode. Given the "shuttle-like" motion of the CM mode, one 
could expect that both sites populate at the same time when they approach together the 
left electrode and then depopulate when they approach the right electrode. On the other 
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FIG. 10: Hybridization scheme for ABL mode (a) and for CM mode (b) considering only the two 
first vibrational levels n = and n = 1. The levels on the left and right are without electron- 
vibration coupling, but the right levels are shifted by one quantum of vibration with respect to the 
left levels. The coupled orbitals have a different vibrational state n because they are mixed by the 
electron-phonon interaction. It is clearly seen that in the ABL case, the coupling does not mix the 
molecular orbitals between them (<7 n=0 -<7 n= i, cr* =0 -o"* =1 ), while in the CM case, the orbitals are 
mixed (cr n=0 -a* =1 , cr* =0 -a n =l)- 

hand, the ABL motion is a stretch mode and hence, one site would populate while the other 
one depopulates. 

However, in Fig. [TT] (a) we see that the ABL mode leads to population of both sites at 
the same time, while in presence of the CM mode, Fig. [TT] (b), the population sequence 
is shifted and dependent on time. In this case, the wave packet has been centered about 
the energy e a and its energy span is smaller than the difference between e a * and e a . The 
reason for this behavior is that one should think in terms of molecular orbitals rather than 
site orbitals. In this case, the ABL motion does not couple the a* and a orbitals, hence, 
the electron populates one of the molecular orbitals, a*, at one time and interacts with the 
vibration in the same way as in the case of the single site. Instead of a single site, we have 
a single level, otherwise there is no physical difference. Let us approximate the wave packet 
by the evolution on the two sites: 

W)ablI 2 ^I^I 2 
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FIG. 11: Population of the ABL mode (a) and of the CM mode (b) on the two different sites 
that are connected to the vibration as a function of time (atomic units). The ABL mode shows 
the same population for the two sites as a function of time, while the CM mode presents a time 
dependent sequential population. 

without any time dependence (except the one due to the coupling to the electrodes that 
we have neglected in this simplified discussion) because ip a is basically an eigenstate of the 
Hamiltonian. This explains Fig. [TT] (a). 

In the case of the CM mode, we saw above that a single-site analogy cannot be applied 
and we have the two molecular orbitals involved in the wave packet propagation. Let us 
again approximate the wave packet by the evolution on the two sites: 

W)CMl 2 « hM 2 + |VV| 2 + 2Re{r^,e i ^-^ t ) 
where a clear time dependence subsists. In good agreement with Fig. [TT] (b). 

B. Two-sites and two vibrations 

When more than one phonon per mode is needed to be converged with respect to the 
electron-vibration coupling, the inelastic transmission cannot be separated in additive con- 
tributions from different modes and all of them have to be considered at the same time. 
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In the above case, the truncated Hamiltonian reads now: 
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Here the subindex a refers to the ABL mode and c to the CM one. At intermediate 
coupling, we need some 3 phonons to converge the transmission. Hence, the Hamiltonian 
matrix in this basis set is composed of 9 x 9 blocks, each block contains the infinite number 
of electron sites, that we deal with as above. 

The inclusion of several modes changes the transmission from the superposition of peaks 
coming from each mode because the spectral weights change and the peaks shift. This is 
seen in Fig. El There the peaks of the full transmission are shifted with respect to the sum 
of peaks from the CM and ABL modes. This can be understood in terms of the different 
Hilbert space that is considered as new modes are included. Indeed, since the convergency 
in phonons, n, means the shift of peaks, the transmission also needs to be converged with 
respect to modes. In terms of Green's function this is easily understandable because there 
are terms in the self-energy stemming from all modes. 

We can improve our previous hybridization diagram to describe the peak structure, by 
including all modes. This type of correlation diagram has been termed "progression of 
progressions" and used in the literature to explain the vibronic sequence of C 60 and naph- 
talocyanine molecules^^i. One word of caution is important: the knowledge of the actual 
coupling matrices is needed in order to perform the correct hybridization of orbitals. The 
hybridization scheme can be obtained by symmetry arguments alone, however the separation 
and strength of the vibronic peaks need a quantitative evaluation of the matrix elements, 
which is increasingly difficult with the number of involved modes. 

VI. THE MEANING OF DIPS AND PEAKS IN THE D 2 I/DV 2 

The increase or decrease of conductance over the vibrational threshold has been used to 
determine the occurrence of vibrational excitation during electron flow through molecular 
electronic states^ 1 ^. 
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A simple explanation for the increase of conductance in the tunneling regime was already 
advanced in IETS of metal-insulator-metal junctions^. The conductance increases because 
new conduction channels become available above a vibrational threshold. This interpretation 
is correct when the vibrational side bands of the electron transmission lie beyond hQ of the 
electrode's Fermi energy. As we have already seen, opening a new channel means that we 
have to add the n = 1 contribution to the n = contribution of the transmission. Hence, the 
d 2 I/dV 2 will present a positive peak at positive voltage, and negative at negative voltage^. 

However, we can also have decreases of conductance. This is experimentally found in 
tunneling^ and in contact regimes^. In the present theory, one always adds a positive con- 
tribution to the transmission above the vibrational threshold, but the change in conductance 
is given by the slope of the transmission, and this can have a rapid variation in the presence 
of a sideband. In the case where the sideband is in the tail of the main peak, and the width 
of the main peak is in the order of the vibrational frequency, then the slope will change 
from a smoothly varying one to the faster varying sideband slope, giving rise to a negative 
d 2 I/dV 2 at positive voltage. Let us give two examples of the above cases. 

A. Tunneling regime 

In the case of the excitation of the CO frustrated rotation mode^, the 2tt* and the 5a 
molecular orbitals are coupled via the electron-vibration coupling. This is easily seen by 
symmetry arguments^, since the mode is antisymmetric with respect to the planes contain- 
ing the molecular axis and the non-zero matrix elements will couple a symmetric orbital 
(5a) with an antisymmetric one (the 27r*)^. Assuming weak coupling, we can neglect the 
effect of all other modes and estimate the change in conductance by using the two-site model 
coupled to a vibration that we presented above. Figure [12 shows the result of the conduc- 
tance, Eq. ( ITOl) . and in the corresponding d 2 I/dV 2 . We see that the Fermi energy plus Q for 
the frustrated rotation, gives the threshold where the transmission for n — 1 enters, giving 
a clear discontinuity in the conductance (positive contribution) and hence a positive peak 
in the d 2 I/dV 2 . This is the case in metal-insulator- metal junctions^. The two-level model 
leads to different inelastic efficiencies at positive and negative voltage bias. 

However, the case of IETS of O2 on Ag (110) 62 is different. Experimentally, dips are found 
instead of positive peaks in the d 2 I/dV 2 . Let us consider the excitation of the internal stretch 
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FIG. 12: Conductance for a model system of two levels of CO on Cu(100). In the inset the 
d 2 I/dV 2 is depicted. The asymmetry in the size of the change in conductance is a consequence of 
the two-level model. The variation in the conductance away from the vibrational threshold is due 
to the energy dependence of the two-level density of states. As can be seen in the (Pl/dV 2 , the 
variation is small with respect to the inelastic change in conductance. 

mode. By computing the PDOS on molecular orbitals^, we know that the n g perpendicular 
to the surface is close to the Fermi energy, and that tunneling takes place through it^. The 
mode has the same symmetry as the molecular orbital and the diagonal matrix element will 
be different from zero. We are hence in a case that can be approximated by the single-site 
model. Figure [13] shows the result. We have located the n g at the Fermi energy as the 
PDOS on molecular orbitals^ seems to indicate. The orbital width is taken as 0.1 eV and 
the frequency energy is 0.08 eV. We see that the transmission function is an asymmetric 
peak because the sideband is inserted in the tail of the main peak. We also include the 
contribution of n = 1 above the Fermi energy plus the frequency in order to calculate the 
d 2 I/dV 2 , as for the CO case while we prefer to plot the complete one-electron transmission 
because in the present case the initial level occupation matters and we cannot treat it in this 
one-electron approach. Instead, we realize that the see the change of slope due to the side 
bands in the n = and n — 1 contributions. The d 2 I/dV 2 gives a dip corresponding to the 
change of slope of the n = contribution. Hence, the elastic component of the transmission 
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FIG. 13: Tranmission for a model system of two levels of O2 on Ag (110). In the inset the d 2 I/dV 2 
is depicted. The decrease in the d 2 I/dV 2 originates in the change of slope of the transmission due 
to the vibrational side bands. These side bands belong to all vibrational channels, in particular to 
the elastic one. Hence, the contribution of the elastic channel alone gives a decrease in conductance 
because the density of states is rapidly changing due to the electron-vibration coupling. 

contains the information of the IETS. The decrease of conductance in this case is due to the 
rapid change of transmission already in the elastic channel. It is then a decrease due to the 
variation of the vibronic density of states. 

From these models we conclude that negative peaks in the d 2 I/dV 2 are due to sidebands 
in the elastic transmission. This statement is equivalent to the one found in the pioneering 
work by Davis 67 . Davis realized that a decrease of conductance could be found in some 
especial circumstances. Namely, that the single-site resonance was at the Fermi energy and 
that the resonance width was of the order of the vibrational frequency. These are the same 
conditions as we find. The interpretation by Davis is that virtual phonons were emitted and 
absorbed giving rise to an interference pattern. This interpretation is based on perturbation 
theory, and it just accounts for the vibration's effect on the elastic channel. In our terms, it 
is just the vibrational sideband of the elastic channel. The only cases in which a vibrational 
sideband can yield a decrease in conductance is when the main peak of the transmission is 
at the Fermi energy- 6 ^ (this is half-filling, we discuss below that we cannot have half-filling 
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in our models), and the electronic widths are of the order than the vibrational frequency, 
so as to enlargened and distort the main peak by the sideband. In the case of half-filling, 
the sidebands are symmetric with respect to the Fermi energy^ giving rise to a peak at 
negative voltage and to a dip at positive one in the d 2 I/dV 2 . The voltage of the dip does 
not exactly correspond to the vibrational frequency, since it correspond to the largest slope 
of the vibrational side band, not its maximum. Hence, these results suggest that the voltages 
at the dips are not direct measurements of the vibrational frequencies, while in the case of 
the voltages corresponding to peaks in the d 2 I/dV 2 , are direct readings of the vibrational 
frequencies. 

There is certain confusion in the literature because perturbation theory is currently used. 
In this case the current in the absence of electron-phonon coupling is confused with the 
elastic current. As can be seen in Fig. [3] (a) in the absence of electron-phonon coupling one 
gets a single lorentzian peak (dashed lines), however the elastic contribution is the n = 
curve of Fig. [3] (b). Hence, it is wrong to identify the elastic current with the one without 
electron-vibration coupling. 

B. Contact regime 

Paulsson and co-workers^ have shown that in the case of symmetric contact to the 
electrodes one can continuously pass from a peak to a dip in the d 2 I/dV 2 by increasing 
the transmission probability. At transmission 1/2 the threshold between both behaviors is 
found. In the case of contact, a large density of states is pinned at the Fermi energy. Hence, 
it is T the parameter that controls the coupling to the leads and, therefore, the passage from 
peaks to dips in the change of conductance at the vibration threshold. Recent experimental 



evidence has been reported in 
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Figure [TH shows the behavior of the electronic wave packet after collision with a 4-site 
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The chosen mode 



chain in the contact regime. The parameters are those from Ref. 
is the ABL, that corresponds to the one giving the largest change in conductance in gold 
monoatomic chains 12 . The wave packet is mainly transmitted in the elastic channel as is 
to be expected from the contact regime, however the vibrationally excited wave packets are 



mainly reflected in agreement with the interpretation of Ref. 



12. However, the behavior 



is more complex: electronic wave packets in channels with odd n are reflected, even n are 
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transmitted. In this particular calculation, the reflection of the n = 1 component is due to 
finite-band effects. As we discussed above, a finite band leads to the closing of transmission 
channels when inelastic effets are operative, increasing the electronic wave packet reflection. 
As the number of vibrational excitations increases, the closing of channels corresponds to 
each excitation process, implying an increase of the reflection in each case. As a consequence, 
even excitations mean even number of reflections, leading to an increase in transmission. 
Indeed, we have found this behavior for large transmission even in the single-site case. 

Some of the calculations of Ref. |69j have been performed in the wide-band approximation, 
and hence the closing of channels due to finite-band effects is not present. Hence, we cannot 
conclude that the decrease in conductance leads to an increase in electron reflection such as 
the one found in this section. Instead, we think their calculations correspond to the regime 
where rapid changes in the vibronic density of states are found which correspond to the 
precedent section. 

VII. MANY-ELECTRON AND ELECTRON-ELECTRON EFECTS IN THE 
PRESENCE OF VIBRATIONS 

Our calculations are exact regarding the vibrational dynamics during electron transport. 
The price to pay is the neglect of multielectron effects. The description of multielectron 
effects together with vibrational excitation necessarily implies approximations. A recent 
treatment of both effects is the one by Galperin et al.— . Despite approximations, their 
treatment is correctly evaluating the vibrational dynamics. Indeed, we can reproduce their 
results, and the main difference between both calculations is an absorption shoulder appear- 
ing below the first peak in the multielectron case. According to Ref. 



20 this is due to the 



propagation of holes coexisting with electron transport. This is absent in our one-electron 
treatment. 

Even in the absence of a direct electron-electron interaction in the Frohlich-Holstein model 
used here, Eq. ©, the electron- vibration coupling can lead to effective electron-electron 
interactions, see for example Ref. 



30i . This leads to electron-electron correlations and has 
extensively been studied in the literature in the so-called bi-polaron problem^. This is 
absent in our one-electron approach, but could in principle be treated by using two-electron 
propagations. 
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FIG. 14: (Color online) Electronic wave packets for different vibrational channels at a time instance 
after the elastic wave packet has passed the vibrating chain centered at site 600. Four electronic sites 
are connected to an ABL vibration. This 4-site chain is strongly coupled to the leads, representing 
the contact regime. Hence, the elastic wave packet is basically unperturbed by the chain, but the 
vibrationally excited channels have wave packets that are reflected for odd n, and transmitted for 



even n. 



Perhaps, more important for the description of transport in metallic systems is the absence 
of the initial electron occupation in the present approach. This is easily solved by usual 
non-equilibrium Green's functions approaches^ at the price of simplifying the treatment on 
the phononic description. Indeed, while the present approach treats the nuclear-coordinate 
description exactly within the one-electron problem and the harmonic approximation, NEGF 
approaches need to rely on approximate treatments^^ 1 ^. This situation can be solved 
by using a time-dependent Hartree-Fock approach to generate a many-body wave packet, 
extending the present treatment to a multielectron case. Some encouraging results in time- 
dependent Hartree-Fock and beyond techniques have been presented in Refs. 
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VIII. SUMMARY AND CONCLUSIONS 

Time-dependent wave packet propagations for the study of inelastic effects in electron 
transport can be very interesting because of the physical picture they permit us to de- 
velop as well as the good size-scaling properties that they can have. Indeed, in the case of 
sparse Hamiltonians, the scaling is basically linear with the system size, when Lanczos-like 
propagation methods are used. 

In the present work we have analyzed the electron-vibration problem when a single elec- 
tronic site is connected to two electrodes. We have used a time-dependent description to 
understand the vibrational excitation sequence, the electronic phase shift and the interfer- 
ence patterns. The time-dependent calculation has given us access to the same quantities 
we can calculate using an energy-resolved theory, but with the time-dependent insight. 

We have also studied the case of two sites and two vibrational modes, and we have 
realized that the two-site problem can be simplified to the single-site one when molecular 
orbitals can be used. The electron-vibration coupling will determine when molecular orbitals 
are meaningful. This is seen by changing the electron-vibration coupling matrix to the 
molecular-orbital basis set. Pursuing this idea, we have developed a simple hybridisation 
scheme that permits us to understand the vibronic structure of an electron transmission 
function in terms of the symmetry of the electron- vibration couplings. We have shown, that 
for medium and strong electron-vibration coupling all modes need to be considered. 

We have applied our simple model to get some insight in the case of the vibrational 
excitation of CO on Cu(100) 63 and O2 on Ag(llO) 62 . Due to our one-electron treatment, the 
calculations are missing fundamental ingredients, but they permit us to explain the IETS 
signals in terms of vibrational channels. In the case of O2, the calculation show that the 
elastic channel leads to a decrease in the d 2 I/dV 2 . Contrary to perturbation theory, the 
elastic channel is not the electronic structure in the absence of electron-vibration coupling, 
but the n = transmission function. Hence the elastic channel contains all of the information 
about the electron-vibration coupling. In particular, the elastic contribution will contain 
vibrational satellites. These satellites or side bands have a strong energy dependence that 
naturally leads to dips in the d 2 I/dV 2 . Indeed this means that the density of states that 
should be considered is the vibronic one, this density of states contain rapid variations with 
the energy that lead to decreases of d 2 I/dV 2 . In perturbation theory, one needs to consider 
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the effect of the vibration in the elastic channel to retrieve the effect of the vibration in 
the electronic structure and hence to take into account the vibronic density of states in the 
conductance. 

We have also explored the case of contact with the two electrodes. In this case, it 
has been shown that the conductance should drop at the vibrational thresholds because 
the vibrations backscatter the impinging electrons. Our calculation show that for the first 
excited (and generally odd vibration channel, n) the electronic wave packet is reflected at 
the vibration. However for even n the wave packet is transmitted. This behavior is due 
to the finite electronic band of the present calculation, that leads to the closing of inelastic 
channels increasing reflection. An even number of consecutive reflections lead to increased 
transmission. 

New developments in time-dependent Hartree-Fock lead us to think that a full time- 
dependent treatment of the physics usually explored with non-equilibrium Green's functions 
will be soon available permitting us to develop a different point of view on inelastic effects 
in electron transport on the atomic scale. 
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